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Abstract 

Solutions of vacuum Einstein's field equations describing uniformly 
accelerated particles or black holes belong to the class of boost-rotation 
symmetric spacetimes. They are the only explicit solutions known which 
represent moving finite objects. Their Newtonian limit is analyzed us- 
ing the Ehlers frame theory. Generic spacetimes with axial and boost 
symmetries are first studied from the Newtonian perspective. The re- 
sults are then illustrated by specific examples such as C-metric, Bonnor- 
Swaminarayan solutions, self-accelerating "dipole particles" , and general- 
ized boost-rotation symmetric solutions describing freely falling particles 
in an external field. In contrast to some previous discussions, our results 
are physically plausible in the sense that the Newtonian limit corresponds 
to the fields of classical point masses accelerated uniformly in classical me- 
chanics. This corroborates the physical significance of the boost-rotation 
symmetric spacetimes. 

1 Introduction 

There is only one class of explicitly known exact solutions of vacuum Einstein's 
field equations which represent moving finite sources: it describes "uniformly 
accelerated particles or black holes". In a curved spacetime the "uniform ac- 
celeration" is understood with respect to a fictitious flat background. However, 
these solutions can be characterized fully geometrically since, in addition to the 
Killing vector associated with axial symmetry, they admit a second Killing vec- 
tor which becomes a boost in the flat-space limit. Therefore, in the name 
"boost-rotation symmetric spacetimes" was coined, and this has beeen used in 
the literature thereafter - among others in the second edition of the "Exact - 
Solutions book" [35] , where some of the solutions with boost-rotation symmetry 
are briefly described in Section 17.2. 

The first solutions of this type were obtained by Bonnor and Swaminayaran 
[Tl] and Israel and Khan [55] in 1964. In 1968 Bicak [1] shown that these so- 
lutions are radiative. Together with other boost-rotation symmetric solutions 
discovered later, they are the only explicit examples of spacetimes with grav- 
itational radiation described by a non-vanishing Bondi's news function, they 
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exhibit peeling-ofF property and admit asymptotically flat null infinity J' sat- 
isfying Penrose's requirements, at least locally [10], [8]. In fact, this situation 
may not change soon since it appears extremely complicated to search for such 
explicit radiative solutions with only one symmetry, and among those with two 
symmetries, the boost-rotation symmetry plays a unique role: Theorems can 
be proven that in axially symmetric, locally asymptotically fiat spacetimes (so 
that exists though not necessarily globally), the only additional symmetry 
that does not exclude radiation is the boost symmetrj0. 

In Minkowski spacetime the boost Killing vector associated with the boost 
along the z-axis has the form 

C,_, = z_-,ct-. (1) 

Here we have to keep the speed of light c because of our interest in the Newto- 
nian limit. The orbits of symmetry to which the Killing vector is tangent are 
hyperbolas — c^t^ = B = constant, x, y ~ constant. Orbits with B > 
are timelike. They can represent worldlines of uniformly accelerated particles 
in special relativity. If a point charged (electric or scalar) particle moves along 
the z-axis with a uniform acceleration and is momentarily at t = at rest at 
z = VB^ its worldline is the hyperbola x = y = 0, z = \/c^t^ + B, (? j^fB is its 
acceleration; its (scalar or electromagnetic) field has boost-rotation symmetry. 
The field is analytic everywhere outside the source if there are two particles lo- 
cated symmetrically with respect to the plane z = and accelerated in opposite 
directions, and both retarded and advanced effects are admitted. 

Such particles are shown in space and spacetime diagram in Fig. [T] From 
the figure it is evident that the two particles move independently of each other 
since their worldlines are separated by the null hypersurfaces z = ± ct. Following 
[TT| , we call the null hypersurfaces z = ± the "roof" . The boost Killing vector 
becomes null on these hypersurfaces, it is timclike for (?t^ < z^, i.e., "below the 
roof", but it is spacelike "above the roof" (c^t^ > z^) - cf. Fig. [T] Since 
the particles' worldlines approach null lines both in the past and future, it is 
intuitively clear that they "start" at the past null infinity and "end up" at 
the future null infinity J^^. 

These features of the boost-rotation symmetric fields in Minkowski space 
emerge without a change also in the curved spacetimes with boost-rotation 
symmetry in general relativity. From a unified point of view these spacetimes 
were defined and treated geometrically in [11] (see also reviews [5], [31]). We 
refer to this detailed work for rigorous definitions and theorems, including the 
global analysis and investigation of the properties of null infinity .y. In gen- 
eral relativity, the "causes" of the motion are incorporated in the theory. In a 
general case of the boost-rotation symmetric solutions there occur conical sin- 
gularities of the metric distributed along the z-axis (see Fig. |T]). They can 
be considered as "strings" (or "rods") which cause the particles to accelerate. 



^See 1101 for the precise formulation and the proof of the theorem, and [5] for the general- 
ization to the electrovacuum spacetimes with Killing vectors which need not be hypersurface 
orthogonal. 
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Figure 1: Space (a) and spacetime (b) diagram with two uniformly accelerated 
particles. The "strings" are depicted in (a) by zigzag lines, their history is in 
the shaded region (b). 

They appear also at J^, so that some parts of its generators are missing; how- 
ever, the distribution of nodes/strings can always be arranged in such a way 
that admits smooth regular sections [11]. In the exceptional cases, when 
is regular except for four points, either the particles are "self-accelerating" due 
to their "inner" structure which has to include a negative mass, or there are 
more particles distributed along z > 0, and symmetrically along z < 0, with the 
signs and the magnitudes of their masses and accelerations properly chosen. An 
infinite number of solutions with self-accelerating particles was constructed in 
[7] ; a pair of chasing particles with positive and negative mass was considered 
approximately in the region > c^t^ first by Bondi [12] with the particles 
being "real extended bodies" ; the corresponding complete exact solution with 
particles represented by (Curzon-type) singularities, is included in |14| and its 
radiative properties are analyzed in [J . 

In the exact boost-rotation symmetric solutions available explicitly the "par- 
ticles" or "sources" can be represented not only by singularities but also by black 
holes. One can construct two or arbitrary even number of oppositely located 
uniformly accelerated black holes connected by a string which pulls them apart. 
Then J" admits regular sections and spatial infinity is regular. Alternatively, 
each black hole is attached to a cosmic string extending to infinity, the region 
"between the holes" is regular. In general, conical singularities extend along the 
whole z-axis but that between the holes has a smaller deficit angle; when the 
conical singularities are interpreted as cosmic strings, the difference in string 
tensions provides the cause of the acceleration. All these cases are described by 
the well-known C-metric; see, e.g. [17] j [13] j [IH]j [IS]j and references therein. 
The distribution of the strings and the sources along the z-axis can be inter- 



preted for the particles represented by singularities and black holes in the same 
way. There also exist "generalized" boost-rotation symmetric spacetimes which 
contain neither cosmic strings nor negative masses. They describe accelerated 
particles in asymptotically "uniform" external fields which can be obtained by 
appropriate limiting procedures from asymptotically flat solutions [Hj, All 
the special cases of the boost-rotation symmetric spacetimes mentioned above 
will be considered in the following. 

Although much understanding of the spacetimes with boost-rotation symme- 
try has been gained, some possible roles of these solutions in general relativity 
have not yet been fully elucidated. Various specific boost-rotation symmetric 
solutions were employed as test beds in numerical relativity based on a null-cone 
formalism [9] , [25] , and on a standard spacelike initial hypersurface [T] . Already 
in [4], where the radiative character of the specific boost-rotations symmetric 
solutions was first discovered and in a "weak-field limit" compared with radi- 
ation from an analogous electromagnetic system, it was suggested that these 
solutions could be used as tests of analytical approximation methods. The first 
natural step in a systematic comparison with post-Newtonian approximation 
methods is to investigate the Newtonian limit. 

Numerous useful, though frequently heuristic, studies of such a limit ap- 
peared since the birth of general relativity. The most geometric ~ and from 
various viewpoints most profound - treatment was initiated independently by 
Cartan [T5] and Friedrichs [23] in 1920's with their formulation of New- 
ton's theory in a generally covariant 4-dimensional spacetime language. This 
so-called Newton-Cartan theory led to the Ehlers frame theory [22], [20], [21] . 
The frame theory is encompassing general relativity and Newton-Cartan theory. 
Our treatment is based on this theory. 

The Newtonian limit of the boost-rotation symmetric spacetimes employ- 
ing the Ehlers frame theory was, in fact, recently investigated by Lazkoz and 
Valiente Kroon [28] . They review the Ehlers frame theory in detail, apply it 
carefully to the general boost-rotation symmetric spacetimes, and illustrate by 
specific examples. However, the results they obtain exhibit surprisingly "un- 
physical features" , using their own words: "The interpretation of particles mov- 
ing in a uniformly accelerated fashion is only valid for early times t « 0." Or 
elsewhere: "Thus the Newtonian limit for early times is a strictly static New- 
tonian potential in which the sources are not moving." Moreover "the boost 
Killing vector field is not inherited" in the Newtonian limit, although one would 
expect that it will just go over to the boosts in the Galilei group. The authors 
of [55] attribute these unphysical features of the Newtonian limits to "problems 
already existing in the general relativistic solutions" . 

In the present work we show that one can construct the Newtonian limits 
with entirely plausible properties. How can such seemingly conflicting conclu- 
sions be understood without finding a mathematical inconsistency in one of the 
treatments? The answer lies in the fact that the construction of a Newtonian 
limit is not a uniquely prescribed procedure. First, one distributes "appropri- 
ately" parameter A = c^^ into both the starting general-relativistic metric and 
matter variables. In this way an one-parameter family of spacetimes is con- 
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structed. The Newtonian limit is then obtained by letting A ^ 0. However, 
there exists no a priori prescription for the construction of this family (or, the 
"allocation" of A's at the beginning). In such a situation the best guidance 
is to anticipate a "reasonable" Newtonian limit and to construct the family of 
spacetimes accordingly. 

Now Lazkoz and Kroon [21] motivated their choice as follows: 

"The natural arena for the discussion of the Newtonian Hmit of boost-rotation symmetric 
spacetimes is the region above the roof, where the spacetime is radiative. The reason for this 
is that as one makes A — > 0, the region below the roof gets squeezed by the roof. That is, the 
region below the roof disappears in the limit, while the roof becomes the set {t = 0}." 

Here, it is, of course, correctly pointed out that the null hypersurfaces of 
the roof degenerate into a spacelike surface in the limit A ^ (like null cones 
get flattened with c oo). However, to stay in the radiative part of the boost- 
rotation symmetric spacetimes means to "loose the particles" ! No wonder that 
the resulting limit has unphysical features. 

In our approach we stay below the roof, i.e., in the region where the boost 
Killing vector is timelike and where the sources occur. In order not to get 
"squeezed" by the roof as A ^ 0, we have to "run away" with the sources. In 
accordance with this idea we parametrize appropriately the particles' worldlines. 
In this way we arrive at the Newtonian limit with all the features expectable 
within the Newtonian gravity. 

In Section [2] we first very briefly summarize the Newtonian limit as it is 
formulated in Ehlers's frame theory [21] - [22]. Since the detailed expositions of 
this theory are not easily accessible, the reader may wish to consult the work 
of Lazkoz and Valiente Kroon [28] where the frame theory is well described, 
including rigorous deflnitions, theorems and references to other works related to 
the theory. We then explain in detail our approach to the Newtonian limit of 
general boost-rotation symmetric spacetimes. The Newtonian limits of several 
typical examples of the solutions with uniformly accelerated sources, in partic- 
ular those mentioned above, are studied in Section [3] Some concluding remarks 
comprise Section 21 We also add an Appendix in which the Newtonian limit 
of the Schwarzschild solution is studied within the whole Schwarzschild-Kruskal 
manifold. There are well-known similarities between the Schwarzschild horizon 
and the acceleration horizon (the roof) , so that the comparison of the Newtonian 
limits is illuminating. 

2 The Newtonian limit of the boost-rotation sym- 
metric spacetimes: The general case 

In contrast to [28j , our construction of the Newtonian limit docs not start from 
the set of observers which fill in the region "above the roof" {c^t^ > z^) where 
no timelike Killing vector exists and no sources occur. Rather, we consider all 
timelike worldlines with fixed canonical coordinates {x, y, z} (see for their 
geometrical meaning), i.e., equivalently, with fixed {p, ip, z} related to them 
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as cylindrical to cartesian coordinates in flat space. Our starting form of the 
metric in the canonical coordinates is 

„ e^'c^ [zdt-tdzf - e'' izAz- cHdtf „^ „ „, , 
— c^r^ 

where ^ = n{p^ , c^t^ — z^), ly = c^i^ — 2^) and the coordinate ranges are 

t, z e (—00, 00), p £ [0, 00), </j e [0, 27r). 

This metric is global - it is valid both above and below the roof (not just 
above it as stated in |28|). As a consequence of vacuum Einstein's equations, 
function fi satisfies the ordinary flat-space wave equation 
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pdp 


(4) 


dz^ 



M = 0. (3) 



c2 dt^ 

The other metric functior0 v is determined by nonlinear equations 

{A + B) v^A = B {Afi^-A - Bp% + 2Afi^A M.b) 

+ {A-B) fi^A - 2Bn,B , (4) 

iA + B)iy^B = A{Bp%- Aii\ + 2Bti^Ai^,B) 

+ {A-B) ^l,B + 2A^i^A , (5) 

where A = p'^, B = z^ — c^t^. The integrabihty condition of Eqs. (g]) and ^ is 
just the wave equation ([3]) which, written also in terms of A and B, reads 

Ap.^AA + Bfl^BB + + M,B = . (6) 

Following Ehlers's frame theory, we consider a family of 4-dimensional man- 
ifolds ^(A) parametrized by the "causality constant" A (= l/c^ for A 7^ 0). 
Each manifold is endowed with a torsion- free connection F^^ ("gravitational 
field") and with temporal and spatial metrics which are given by symmetric 
tensors taf3 and s"^. These metrics satisfy the relation 

tcaS^^' - -A(5^ , (7) 

and are compatible with the connection Fg^. The frame theory requires that 
for A 7^ the metrics arc tap = —Xgap and s"^ = g"'^. 

In the limit A ^ the frame theory goes over to the Ncwton-Cartan the- 
ory - a geometrical formulation of the classical Newton theory (see, e.g., [37] 
and references quoted in [12]). Then there exists (at least locally) a Cartesian 
coordinate system such that the metrics have simple forms 

tafj = diag(l, 0, 0, 0) , = diag(0, 1, 1, 1) , 



and the gravitational field can be expressed by means of the affine connection 
as 

' 137 



F2^ = t,^t,^s"*$5, (8) 



■^In 111! and other references this function is usually called A. Wo use u instead of A because 
A is here reserved for the causality constant of Ehlers's frame theory. 
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where t is the global (scalar) absolute time and <I> is the Newtonian gravitational 
potential. 

We shall now derive the relevant quantities entering the frame theory for the 
boost-rotation symmetric spacetimes. First let us write the worldline of a uni- 
formly accelerated particle in special relativity (see, e.g., for a pedagogical 
exposition) in the form 

This describes the well-known hyperbolic motion when the particle starts at 
infinity with velocity approaching c, arrives at z|t=o = c^/d = ^/^9 where it 
stops, and then continues back to infinity, approaching it asymptotically with 
V —I- c; g is its proper acceleration. 




(e) A = 1 (f) A = 1/2 (g) A = 1/4 (h) A = 1/20 



Figure 2: The roof and the particles worldlines for different values of A in 
the global coordinates t ~ z (figures (a) -(d)) and for the same values of A 
in the shifted coordinates t — where C, = z — l/Xg (figures (e)~(h)). The 
corresponding regions are shaded. 



Applying the limit A ^ 0, we immediately observe that as A ^ the null 
cone ~ and so the roof ~ becomes the hyperplane t = 0, and the particle, 
including its turning point, is pushed away to spatial infinity. Hence, in order 
to obtain a nontrivial limit, we have to "go" to infinity with the particle. We 
do this by introducing a new coordinate C by 

^ = C+T^, (10) 
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i.e., we make a A-dependent shift of z. It is illustrated in Fig. [5] how this shift 
works for various values of A. Since in the boost-rotation symmetric spacetimes 
particles symmetrically located along z < occur, their worldlines are also 
plotted. However, notice that these are "shifted away" as A — > 0. 

Wc now apply this procedure to find the Newtonian limit of the boost- 
rotation symmetric spacetimes. Replacing the factor c by 1/-\/A in the metric 
([2|), the metric tensor in coordinates {i, r, z, Lp\ turns out to be 



2A-i2 





(e^ - e'' 



tz 



Z2A - i2 












z^X-t^ 




e^z^X 



z^X - 




2 - 

-re 



(11) 



its contravariant form is 



( ^e-f'z^X - e-^t^ 

' z^X - <2 


tzie-i' - e-") 



V 



'-X~t^ 











tz (e 



Z'^X - i2 



-''z2A-e-^t2 

Z2A - t2 
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(12) 



Factor A~^ appears also in A^^t2 in the arguments of functions /i and v. In this 
way the A-factors arc "distributed" uniquely in the metric and the Newtonian 
limit is thus fixed. Then, as noted above, the temporal and spatial metrics in 
the frame theory arc simply tap = —Xgap and s"^ = g"^ ■ 

By analyzing all specific cases in Section [3] we shall see that when A ^ the 
functions n and v are approaching constants and all their derivatives fi_a and 
v,a vanish. This can be seen on physical grounds by regarding Eq. ^ with a 
sourc^, i.e. 

□Ai = SirGXa . (13) 

'1 r_i _ i.„ .^-3 



So, n is 0(A) plus possibly a constant. (Since [G] 



•kg 



kg-m 



[A] = m~2 . s2, is indeed dimensionlcss.) However, by adding constants to ^ 
and v in exact solutions one introduces or removes conical singularities 
("strings" or "struts") along the z axis (see, e.g., [H]). These, in generic cases, 
are attached to particles and cause their motion, as indicated in Fig. [TJ If 
they should not represent just conical singularities introduced "ad hoc" , they 
turn out to be proportional to the products mA, the mass and the acceleration 
parameters of the particles (see examples below) . These products are ^ A2 and 
can thus be neglected in the Newtonian limit. Hence, we can assume fi ^ 0(A) 



•^One of the Ricci tensor components turns out to be R^^ ~ 
hand side of the Einstein's equation becomes iT^iptp ~ f gipipT) 



' "p-^ D/x, and the right- 



^ p^e f^XGcr for an ideal fluid. 



8 



and V ^ 0(A) without additional constants. The axis is thus regular in the 
Newtonian limit, except for places where the particles occur. 

Before performing the limit A ^ it is crucial that we make the substitution 
(fTU|) . Then wc find that, as A ^ 0, the metrics tap and s"'' approach 

tap = diag(e^ 0, 0, 0) , s"'' = - diag(0, e'^ e'^ / p") , (14) 

where jl = Xmix^Q ii{p^ , (C + ^)^ — t^/X) and analogically v. These limits 
indeed give p ^ v = Q, as we discussed above, but we leave metrics in the form 
(fT4|) to show explicitly how substitution (fTO)l alters p and v after the limiting 
process. 

Starting out from (lllll and p^ . we can calculate the affine connection and 
so find 

1 y".a' 



T% = lim 



2 A 



{p, z}, (15) 



(there are also non- vanishing components and T^^ arising due to the cylin- 
drical coordinates - these are not needed here). Notice that at this point the 
metric function v no longer enters the Newtonian results. 

As a consequence of Eq. ([TS|) wc identify the Newtonian gravitational po- 
tential as 



M(p^(c+i)'-^VA) 



$ = lim — ^ ^ . (16) 

(No "minus" sign arises here because of our signature.) This is consistent with 
the affine connection ([5]) in the sense that we obtain ([TS]) from with t and 
s given by ([14]). 

In fiat spacetime with Minkowski metric (1/A, —1, —1, —1), two symmetri- 
cally located point particles move along the worldlines = {t{T), 0, 0, Zeir)) 
with ^ ^ 

Zeir) = e — cosh ^VAgr^ , t{T) = —j=- sinh {VXgr^ , (17) 

where e = ±1 denotes the particular hyperbola. Their 4- velocities are normal- 
ized to 1/A with respect to Minkowski metric. The matter densities are thus 
given by (Te = / (5'-^-'(a;'' — x^) dr. After a short calculation wc obtain the total 
density 

a = (T+l +(T_l=— o z — 



which, after substitution (jlOp and limit A ^ 0, yields a = 5{( — ^gt'^), i.e., we 
arrive at the expected Newtonian result for the source given by one uniformly 
accelerated particle. 

Before turning to the Newtonian limit of the field, we should prove that 
this form of motion is not merely a peculiar coordinate artefact. We shall show 
that the motion takes place with respect to privileged ("rigid") observers. We 
started with the "observer field" u^^ = (l/y^ioo, 0, 0, 0) which is normalized 
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to 1 with respect to the temporal metric ta/3 (or to 1/A with respect to ga/s)- 
Following the Theorem 2.2 and Corollary 2.3 of [28], the important information 
is contained in the tensor 

= hm \ ur. J , or F = lim \ du' , 

where u' = taf}U°'dx^ , the bullet indicates that indices arc lowered by temporal 
metric tap- The observer field has the limit = (1, 0, 0, 0) which remains nor- 
malized with respect to the temporal metric. This represents the field measured 
by static (test) observers fiowing only in the absolute time t. They perceive a 
gravitational force F = — V<i> caused by potential $ which depends on z — ^gt^ 
and thus represents the potential of a moving source. 

We shall now demonstrate that in the limit A ^ the field equation p3p 
implies the standard Poisson equation for the potential $. Introducing new 
variables a = A and b ^ B —jr^ = — t^ /\ — we find Eq. with a 

source a = a{a, b) to becomcl 

a^Ji.aa + tJ.,a+ {b+ -gTj2^ ^^^.bb + lJ-,b = SttGAct , (18) 

where d^A — d,a and d,B = d,b because the variables differ just by shift pa- 
using d_b = A9 At, and going back to z and t, the last equation reads 

aH,aa + Ai,a + A (z^A - t" ) 



d 



+ X—, — ^ = 27rGAcr. (19) 

Taking now the substitution (|10[) into account, we find the limiting values of 
the factors entering this equation as follows: limA^o A {z^X — t^) = and 
\imx-,o{z^X - ~ -^) I (C- \gt^)- Dividing (HH) by A, setting, $ = 
lim^^o /i/2A following (|16p . and putting S = limA^o we find Eq. (|19p in 
the limit A ^ to go over to 

2 (a$,,a + $,„) + \ ^ ^"^ . = 2^GE . (20) 

After introducing finally C = C ^ \9^^ writing back a = p^, we recover 
standard Poisson's equation: 



*The fiat-space d'Alembert operator (|3]l is in terms of coordinates {A, B) given by □ 
4(dA+AdAA + du+Bd,BB)- 
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where $ = $ (p, C) and E = E (p, () . 

These results are in strong contrast with those obtained in [28]. There the 
Newtonian "sources" - and thus also potentials - depend only on p^, t"^ but not 
on z; hence the "sources" are cyhndrical as, in fact, admitted in [55|. 



3 Examples 

In this section we wish to illustrate previous results on some typical classes 
of the boost-rotation symmetric spacetimes. We give the Newtonian limits for 
two accelerated "monopolc particles" , for Bondi's chasing freely falling particle^ 
with a negative and a positive mass, for accelerated "dipole particles", for freely 
falling particles in an external field, and for the C-metric. 



3.1 Accelerated "monopole particles" 

These are described by the Bonnor-Swaminarayan solution [14) given by metric 
© with 



2AGmi 2XGm2 
AiRi 



/ m 



A2R2 
2 



AXGmiAi + A\Gm2A2 



(21) 



^1 , J!^^^2/2 ,2/^^ _2miTO2 



R2 



-I-8TO1TO2 



AlAl (i?i - R2f 



-2XG 



(^2 - AlY R1R2 

niiAi , 7712^2 



i?2 



+ z'- e/x) , 



where 




Aj 



Aj 



(22) 



(23) 



In functions ([21]) - ([23]) we already inserted A and G so as they have the correct 
form in standard units. In addition, we have to make substitutiorlfl Ai = Xgi, 
where accelerations gi have usual dimension m • s"^. Since now dimensions are 



[.9.] = 



■ s 2 and [G] = m'^ • s ^ • kg 



the 



[Ri\ = m^, [mi] = kg, [A] = m ^ 
function fi is dimensionlcss. 

This metric represents, in fact, two pairs of particles. Only two symmetrl 
cally located particles accelerated uniformly in opposite directions with acceler- 
ation Ai are obtained by putting 7712 ~ 0. 



■"For simplicity and as is common in literature, we are using the term "particle" through- 
out this paper although for the static Curzon-Chazy particles the "inner" structure of the 
singularity is quite complicated 1341 (this will not change by accelerating them). Thus they 
are not point-like particles in an ordinary sense. 

^Thesc yl's are acceleration parameters of the solution, not the variables as in equations 
^ — JSjl . It should be clear from the context what is meant by A in a particular equation; in 
this section A denotes only accelerations. 
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Regarding now the general result (jl6p for the Newtonian potential and using 
for its evaluation ((2T|) and (|22p . we find 



$=lim^ = -^^^L 

A^o 2A / o . 



v/p^ + (C-i5it^)' 



This is precisely the potential due to a single point particle with mass mi that 
moves with constant acceleration gi. 

In the limit A ^ 0, function ly and its derivatives go to zero as they should; 
in contrast to /.j, function ly does not enter the Newtonian limit (cf. text below 
Eq. (fTS]) ). Recall also (cf., e.g., Fig. [2]) that our limiting procedure removes the 
symmetrically located ("left") particle to C ^ — oo. 

The case of two pairs of particles is more complicated. The particles have 
different accelerations and the substitution ([TO)) does not stop the second particle 
to run away to infinity as A goes to zero. However, we can keep the particles 
of the "right" pair at constant mutual distance in the hypcrplane t = by 
introducing a length scale L. This can be done by setting 

91 

In this case, we find that function /i in Eq. ([2T|) with both mi and m2 non- 
vanishing, substituted into Eq. (|16p . implies the Newtonian potential to be 

^ ^"^^ ^ ^"^^ (25) 

This represents the potential of two particles with masses mi and m2 undergoing 
the same acceleration 51, and moving on parabolae differing only by the shift 
L. 

There exist an interesting situation in the case of two pairs of particles, 
considered first by Bondi [T^], when in the full relativistic solution contained 
in Eqs. (PTjl - (QH) no singularities along the z axis are present. This can be 
achieved by choosing masses to be 



"^1 = ^,0 c, 9 , m2 



The whole sequence of spacetimes ^(A), (?q/3(A) has the z axis regular. In the 
Newtonian limit we again obtain the potential (|25p but now, as a consequence of 
Eq. (|24p . the masses are mi = L^gi/G and m2 = —L^gi/G, i.e., their magnitude 
is the same, only signs differ. The accelerations are again same. This is the exact 
solution of the Newton's equations of motion for a pair of particles with opposite 
masses which accelerate themselves uniformly in a straight line. 
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3.2 Freely falling particle in an external field 

Rather surprisingly, by a similar method we are able to find the meaningful 
Newtonian limit even in the case of generalized Bonnor-Swaminarayan solu- 
tions which are not asymptotically flat. These solutions were constructed in 
[B] by removing the "outer" particles described by the solution (^J) - (HH) to 
infinity, while increasing at the same time their mass, so that there remain only 
the "inner" - causally independent - two particles, each of them falling in an 
external field. There are no singularities along the z-axis, except for the places 
where these particles occur. Removing the "outer" particle as in JSj/, but starting 
from Eqs. 1121]) - \2}3fl in which constants G and A are inserted in appropriate 
places, we find the generalized Bonnor-Swaminarayan solution to he given by 
the following functions fi and v: 



Making again substitution A — !■ Ag, wc find that our general result (|16p implies 
the Newtonian potential in the form 



The potential does not now decay at infinity, however it can, indeed, be inter- 
preted as the potential of a point particle with mass m moving in a homogeneous 
external gravitational field of strength g (purely time-dependent last term does 
not, of course, contribute to the field strength). Hence, a solution of Einstein's 
field equations representing a moving object in an external field leads, in the 
Newtonian limit, to a simple "intuitive" solution of Newton's theory. 

3.3 Self-accelerating "dipole particles" 

This is the solution obtained from Bonnor-Swaminarayan solutions by consid- 
ering two particles (with masses of opposite sign) of each pair with conical sin- 
gularities located only between the particles. A limiting procedure then brings 
the particles in each pair together and, simultaneously, the magnitudes of the 
masses are increased [7]. (A similar procedure produces standard dipoles in 
electromagnctism.) The resulting spacetime contains no conical singularities, 
just two causally independent self-accelerating dipole particles. The metric is 




$ = lim = — 
A^o 2A 
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given by ([2]), in which 



S>XG 



R 4A2i?3 

2 



P 



(28) 



((p^+z^-tVA)^ 



1 

l4 



4i?3 



+ z2-tVA) (^p2-z2 + fVA + -^ 



(29) 



A is the magnitude of the acceleration of the symmetrically located dipoles, ^ 
is their constant dipole momenlj^ with dimension \^\ = kg • m. To obtain the 
Newtonian limit, we again put A ^ \g and, after using the A-dependent shift 
(fTO|) of the z coordinate as in Eq. (fT6|) . we obtain the resulting limit in the form 



$ = lim -V = 
A^o 2A 



3/2 ■ 



This is the Newtonian potential of a uniformly accelerated dipole. 



3.4 The C-metric 

Let us turn to the last and the most elaborate example - the C-metric. The 
C-metric can be interpreted as two uniformly accelerated black holes with the 
accelerations caused by cosmic strings. It is an example of (perhaps most fre- 
quently used) boost-rotation symmetric solution (see, e.g., [3], for a number of 
references on the C-metric). In the canonical boost-rotation coordinates, it is 
described by "accelerated" rods, representing the event horizons of black holes. 
In metric ([2]), functions v are now given by [13j 



R\ + R2 



1 _ 1 

In )'"^ '^^( , (30) 



1 _ 1 

2Ai 2Ai 



In , „ „ h In ' 



4.R1R2 \Ri + Zi 

{p' + z'- tyX) R, + 2p' (z^ - tVA) + (z^ ^ f-/X p') Z, 

{p^ + z^-t^/x)R2 + 2p^{z^-tyx) + {z^-tyx-p^)Z2 ' ^ ' 



^In the notation of the present paper, the dipole moment is defined hy = 2{rn\/A\ + 
"12/^2). Assuring the correct position of the struts, the dipole moment reads S) = 
2mi (A| — A^) /AiA|. In the limit A2 — » Ai the mass parameter mi has to be rescaled 
to keep the dipole moment ^ constant. 
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where Ri are defined by ([^5]) and (i = 1, 2) read 



The rods in tlic hyperplane t ~ are segments on the z-axis given by intervals 
(iMi, I/A2) and (-IM1, -IM2), ^1 > A2. 

Since the Newtonian limit of a Schwarzschild black hole is a point mass, 
we anticipate the Newtonian limit of the C-metric to be an accelerated point 
particle. The rod should thus shrink to a point. This is achieved by putting 

A1.2 = , ^ . (32) 

/ 1 T 2mG 

Fa^ + g 

This is more complicated substitution for the acceleration parameter than in 
the previous cases. However, notice that for small mass parameter m we again 
get A ~ gA as before. The parameter m enters the expression for Ai_2 since it 
is related directly to the "length" of the accelerated rods. Recall that in the 
static case in Weyl's coordinates the Schwarzschild black hole is represented by 
a rod of the length 2m where m is the Schwarzschild mass. 

After long but straightforward calculations of the limit p^ . we obtain the 
resulting Newtonian potential. Remarkably, complicated form of the rclativistic 
C-mctric leads to the simple classical potential of uniformly accelerated point 
particle with mass m: 

Gm 



4 Conclusions 

We have constructed physically plausible Newtonian limit of general asymp- 
totically flat boost-rotation symmetric spacetimes. We have also seen that 
a satisfactory Newtonian limit can be obtained in some cases which are not 
asymptotically flat. 

In case of rclativistic spacetimes without conical singularities, when parti- 
cles fall freely under mutual gravitational interaction or in an external field, the 
Newtonian limit produces exact self-consistent solutions of the field (Poisson) 
equation and of the equations of motion in Newton's theory. This satisfying re- 
sult increases significantly the physical meaning of the boost-rotation symmetric 
spacetimes. 

There exist other solutions for which we can apply the above method of 
finding the Newtonian limit and so extend our results. For example, one can 
explicitly write down functions ji and v representing two (or even more) pairs 
of Schwarzschild black holes, Curzon particles, etc. [35]. Also, the solution 
representing a Schwarzschild black hole moving in an external field is available - 
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the "generalized" C- metric [23] , which is the analogue of the generalized Bonnor- 
Swaminarayan solution discussed in Section [3.2l Their Newtonian limits can be 
expected to be of the same form as those obtained above. 

In our procedure of performing the Newtonian limit, we lost the track of 
the three quadrants of the original relativistic spacetime, namely, c^i^ > and 
(?t'' < for z < 0; hence, the particle(s) accelerated in opposite direction 
along the negative part of z-axis "disappeared" in the limit. In Appendix, we 
indicate that an analogous result arises in the Newtonian limit of the complete 
Schwarzschild-Kruskal spacetime. In fact, the region "preserved" during our 
limiting procedure is that in which, in the original spacetime, the boost Killing 
vector field is timelike and the moving sources occur. From the Newtonian 
perspective, this region is far more intuitive than the "radiative" region above 
the roof considered as "the natural arena" for discussion of the Newtonian limit 
in [28] because "the region below the roof gets squeezed by the roof" . 

The existence of physically "correct" Newtonian limits of the boost-rotation 
symmetric spacetimes, in particular of those in which only particles/black holes 
occur without any strings or struts, makes it plausible to assume that the exis- 
tence of complete global relativistic solutions can be proven, at least when the 
fields will be "close" to their Newtonian counterparts. Sources will be repre- 
sented by non-singular extended bodies. Bondi [12] discussed spatially extended 
material "balls" of a positive and negative mass uniformly accelerated by mu- 
tual gravity in the region below the roof, but did not succeed in constructing 
exact solutions. It is conceivable that when elastic bodies are considered, the 
existence of near-Newtonian relativistic solutions can be proven [33j . 

In this journal, Antoci ei al. |2] recently questioned the physical meaning 
of the boost-rotation symmetric spacetimes, arguing that the acceleration hori- 
zons in these solutions (the "roof" in our terminology) are singular. In their 
Conclusion, they also found "problematic" the common interpretation that the 
singularities in these solutions represent the masses exhibiting a uniformly ac- 
celerated motion. Their claim concerning the acceleration horizons was shown 
to be unfounded in a detailed analysis of MacCallum [23 . The interpretation of 
singularities as accelerated sources can be supported in various ways (see, e.g., 
[TT]). An important argument for this interpretation comes from the Newtonian 
limit considered above. 

We are grateful to the Albert Einstein Institute, Golm, for the kind hos- 
pitality, and the Grants No LC06014 ("Center of Theoretical Astrophysics") 
and MSM002 1620860 of the Ministry of Education of the Czech Republic for 
the partial support. J.B. also acknowledges the support of the Alexander von 
Humboldt Foundation and from the Grant GACR 202/06/0041. Last but not 
least, J.B. thanks Jiirgen Ehlcrs and Bcrndt Schmidt for useful discussions. 
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A Appendix: The Newtonian limit of the Schwarzschild- 
Kruskal solution 



We sketch the Newtonian hmit of the Schwarzschild solution in the global 
Kruskal coordinates. Our Newtonian limit of the boost-rotation symmetric 
solutions and the character of Killing horizons show some similarities. The 
Schwarzschild metric in standard coordinates is well suited for the Newtonian 
limit - see |20| : 

d.^ = 4 (l - ^) dt^ + ^-4^ dr^ + . (33) 

Let us remind the standard procedure of introducing global Kruskal coordinates 
[30] . keeping G and A = c~^. Introduce tortoise coordinate r* by 

.*=. + 2.nAGln(^-l), (34) 

and then the null coor dinates U and V by V ~ U 2r* and V + U ^ 2t/y/\. 
Next, coordinates u and v are given by u = exp(— 4^J^xg ) ^' ~ exp(— ^rnXG ^' 
Finally, introduce global time T and coordinate Rhy v ~ u = ^^^q and v + u = 
Jfe- Then 



These coordinates have correct dimensions of length and time. 
The Schwarzschild metric ([55)1 becomes 

ds^ = ^^^^e~^ ("l*^^' + '^^') + ■ 

It is invariant under Lorentz transformation mixing coordinates (i?, T) together 
since i?2 _ T'^jx = m^A^G^ [r/2m\G - 1] exp (r/2mAG) is unchanged. 

The inverse of relations ([55]) - ([35]) cannot be expressed explicitly, so r is to 
be understood as the implicit function of R and T. It is useful to express the 
inversion in terms of Lambert W functiorlf] as 

r = 2mAG 
t = mA'^^^G arcsinh 



* Lambert 'P' function is the solution of the equation x 
about Lambert W function see 1171 . 
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The Schwarzschild metric ([57)) takes the form {w denotes the parameter of 
Lambert W function, i.e., w ^ {R'^ - T^/A) /m^A^G^e = -uv/e): 
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H 1 



e W{w) + lw \ A 



1 



-- dT^ + di?^ + 4?7i^A'G-' + ly dn'' . (40) 



This is wholly expressed in terms of Kruskal coordinates. For all possible values 
of R and T, the value of w lies in the interval [—1/e, oo) where the Lambert W 
function is well-defined. 

Using the relation ^ = ^'^('^)' ^"^^ ^ ~ ~uv/e, we can write 



ds^ 



16 

e 



9w 



A 



dT^ + di?-" + 4m^A^G'' [W{w) + If dfl' . (41) 




(a)A = l (b)A = l/2 (c)A = l/4 (d) A = 1/20 



(e) A = 1 



(f) A = 1/2 



(g) A = 1/4 



(h) A = 1/20 



Figure 3: The Newtonian limit of a Schwarzschild black hole in the Kruskal 
coordinates {T,R) (figures (a) -(d)) and in the original - from the Newtonian 
point of view more physical - Schwarzschild coordinates (t, r) (figures (e) - (h)). 
This should be compared with Fig. [2] in the main text. 



The analogue of the roof gets "squeezed" in the limit A ^ 0, whereas world- 
lines r = const, (i.e., hyperbolas) in the Kruskal coordinates (T, R) run away 
to infinity as can be seen from Eqs. (|35p - ([55)1 . The worldlines of constant r are 
worldlines of uniformly accelerated observers whose acceleration is set precisely 
to balance the attractive force of the central black hole; these worldlines are 
tangent to the Killing vector d/dt. 

The structure of the Killing vector d/dt in the Kruskal coordinates in the 
vicinity of the event horizon resembles the boost vector in the vicinity of accel- 
eration horizon: ^'^ ~ {R^/X, T/^/X) (although topologies differ). 
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